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1.  Introduction 


Advanced  sensors  and  computer  models  for  the  surveillance,  detection,  identification,  and 
tracking  of  sound-emitting  targets  in  combat  are  rapidly  emerging  (1,2).  Most  new  Army 
acoustic  systems  use  unattended  microphone  sensors  to  construct  small  ground-based,  beam¬ 
forming  arrays  to  determine  line-of-bearing  angles,  e.g.,  the  remote  netted  acoustic  detection 
system  (5).  At  the  same  time,  point-to-point  transmission  of  acoustic  signals  is  greatly 
influenced  by  variations  in  the  wind  flow,  temperature  field,  and  turbulence  over  varying 
topography  and  surface  energy  budgets  (4-6).  Thus,  the  Army  is  looking  to  implement  the  best 
possible  acoustic  propagation  models,  particularly  those  that  will  incorporate  the  impacts  of  the 
environment  (e.g.,  the  atmosphere,  turbulence,  and  terrain)  on  acoustic  signatures. 


A  key  element  in  determining  point-to-point  acoustic  transmission  is  modeling  the  variation  of 
the  effective  speed  of  sound  through  the  lower  atmosphere.  The  effective  speed  of  sound  (ceJf), 
which  accounts  for  the  bearing  of  the  sound  source  from  a  given  receiver  and  the  compass 
bearing  of  the  mean  wind,  is  readily  derived  from  measured  or  modeled  estimates  of  air 
temperature,  relative  humidity,  and  wind  velocity.  It  is  useful  to  define  the  effective  speed  of 
sound  (ceff)  from  the  following  expression  given  by  Noble  (7)  and  Osteshev  (5): 


-eff  co 


+  u  cos(0w  —  71  —  Or  )  , 


in  which  ca  =  J— 


^  5  r  =  8314.32  J  mol  lKTi  is  the  universal  gas  constant,  T  is  air 
M 


^  n  _ 

temperature  in  Kelvin,  M  is  molecular  mass,  ys  =  —  is  the  ratio  of  specific  heats,  u  is  the  mean 

cv 

of  the  horizontal  wind,  0W  is  the  bearing  of  the  wind  from  North,  Or  is  the  bearing  of  the  receiver 
from  the  source,  and  u  cos (0W  -  n-  Or)  is  the  component  of  the  sound  speed  along  the  direction 
of  propagation  from  the  source  to  the  receiver.  As  discussed  in  Osteshev  (5),  the  effective  sound 
speed  in  equation  1  is  valid  only  for  nearly  horizontal  propagation  angles  (5).  In  addition,  Wong 


and  Embleton  (9)  have  deduced  the  ratio  of  specific  heats  and  molar  mass 


as  a  function  of 


temperature  and  humidity  in  the  following  form: 


=  0.04833  +  (Rh  -  0.023 )At  ,  (2) 

M 

in  which  Af  =  9.2  x  10~5  +  5.5  x  10~6 t  +  4.25  x  10~7  r,  R,„  is  relative  humidity  (in  the  range  0.0  to 
1.0),  and  t  is  air  temperature  (°C).  Then,  for  applications  in  outdoor  acoustics,  sound  speed 
differences  across  vertical  layers  will  generally  cause  acoustic  waves  to  be  refracted  upward  if 
the  effective  sound  speed  decreases  with  height,  and  they  will  be  refracted  downward  if  sound 
speed  increases  with  height. 
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Monin-Obukhov  (M-O)  similarity  models  (10-14)  have  been  very  well  suited  for  such  acoustic 
applications.  For  example,  the  energy  budget  and  meteorological  models  implemented  in  the 
Acoustic  Battlefield  Aid  (ABFA)  (15,16)  are  relatively  simple  and  quite  robust  for  generating 
low  atmosphere  profiles  of  effective  sound  speed,  especially  for  daytime  convective 
atmospheres.  However,  for  the  nighttime  case,  the  model  in  ABFA  assumes  that  the  net  surface 
radiation  is  constant  with  time,  reduced  only  in  the  presence  of  clouds,  and  is  balanced  entirely  by 
the  surface  heat  flux.  Therefore,  in  this  report,  we  are  interested  in  testing  the  effectiveness  (for 
this  application)  of  an  earlier  quartic  solution  for  the  surface  temperature  (17,18),  i.e.,  an 
alternate  energy  budget  model,  which  explicitly  computes  the  net  (long  wave)  radiative  flux  and 
the  turbulent  fluxes  of  heat  and  moisture  close  to  the  ground.  Then,  the  primary  energy  budget 
model  results  can  be  used  to  estimate  profiles  for  wind  speed,  air  temperature,  and  relative 
humidity,  from  which,  effective  sound  speed  information  can  be  generated.1 


2.  Model  Equations 


2.1  A  Quartic  Solution  for  the  Surface  Temperature 

The  alternate  surface  energy  budget  model  described  in  this  section  is  a  semi-empirical, 
equilibrium  model  (17,18)  that  explicitly  computes  each  energy  balance  component,  except  the 
flux  of  energy  into  and  out  from  the  ground,  which  is  computed  simply  as  a  percentage,  10% 
daytime  and  5%  nighttime,  of  the  computed  net  radiation  (21).  The  model's  primary  input  are 
reference  level  (1  to  2  m)  values  for  wind  speed,  air  temperature,  relative  humidity,  atmospheric 
pressure,  surface  thermal  emissivity,  soil  water  content,  and  cloud  cover  information.  (For 
daytime  calculations,  the  day  of  interest,  longitude,  and  latitude  of  the  site  are  required  to 
determine  the  solar  declination,  hour,  and  zenith  angles.)  The  model's  primary  output  are 
calculations  of  the  net  radiation,  the  effective  surface  temperature,  and  the  soil,  sensible,  and 
latent  heat  fluxes. 

In  accordance  with  Carson  (22),  it  is  assumed  that  the  energy  budget  close  to  the  ground  (as 
shown  in  Figure  1)  can  be  written  as 

Rnet  =H  +  L'E  +  G  ,  (1) 

in  which  Rnet  is  the  net  radiative  flux,  H  is  the  turbulent  sensible  heat  flux,  LE  is  the  latent  heat 
flux  attributable  to  surface  evaporation,  and  G  is  the  flux  of  heat  into  or  out  of  the  soil.  Positive 
values  of  Rnet  and  G  are  chosen  to  represent  fluxes  directed  downward  while  positive  H  and  L’E 
values  represent  fluxes  directed  away  from  the  air-soil  interface. 


'Note  that  Swanson  and  Reichard  (19)  and  Swanson  (20)  have  also  reported  on  the  use  of  M-O  similarity-type  formulations  to 
generate  low  atmosphere  sound  speed  profiles.  However,  they  based  their  estimates  of  the  M-O  scaling  constants  on  data 
retrieved  from  a  ground-based  micrometeorological  instrument  rather  than  from  energy  budget  principles,  as  discussed  in  this 
report.  Nevertheless,  their  goal  was  very  similar,  i.e.,  to  incorporate  similarity-based  estimates  of  effective  sound  speed  within 
probability  of  detection  (sensor  performance)  algorithms  for  unattended  acoustic  sensor  arrays. 


2 


WET  GROUND 

Figure  1 .  The  radiation  and  energy  budget  close  to  the  ground  (illustrated  by 
F.V.  Hansen,  1993). 

The  formulations  in  the  model  used  to  compute  the  incoming  short  wave  solar  radiation  (Rsi)  for 
daytime  cloudless  skies  are  patterned  after  Meyers  and  Dale  (23)  and  are  augmented  as  needed 
with  empirical  results  by  Haurwitz  (24)  to  account  for  cloudy  skies,  i.e., 

Rsl=  I  =  IoTrT  pT\yT  acosZ  ,  (2) 

in  which  I„  is  the  extra-terrestrial  flux  density  at  the  top  of  the  atmosphere  on  a  surface  normal  to 
the  incident  radiation,  Z  is  the  solar  zenith  angle,  and  7)  are  the  transmission  coefficients  for 
Rayleigh  scattering  (R),  absorption  by  permanent  gases  (P),  water  vapor  (W),  and  absorption  and 
scattering  of  aerosols  (A).  The  empirically  derived  transmission  coefficient  for  cloudy  skies  (24) 
is  defined  by  the  ratio  of  insolation  with  partly  or  completely  covered  sky  to  insolation  of 
cloudless  skies. 

Gates  (25)  gives  the  expression  for  downward  long  wave  radiation  (Ru)  for  clear  skies  in  the 
model,  i.e., 

Rl1  =  _  170.9 +  1.195  a  7)4  ,  (3) 

2  —8  —2  —4 

in  which  Rlj,  is  in  watts  per  meter  ,  o-  5.6697  x  10  W  m  “  K  ,  and  T,  the  reference  level  (1  to 
2  m)  temperature  in  Kelvin.  From  Paltridge  and  Platt  (26)  an  addition  to  equation  3  is  made  to 
account  for  clouds,  giving  the  total  expression  as 

Rl1  =  -170.9  +  1. 194  crr,4  +  0.3  scaT^(cc)  ,  (4) 

in  which  ec  is  the  long  wave  flux  emissivity  of  the  cloud  base,  T,  is  the  temperature  of  the  cloud 
base  in  Kelvin,  and  cc  is  the  fractional  cloud  cover.  Having  an  estimate  of  the  cloud  height,  one 
can  approximate  Tc,  assuming  an  average  of  the  dry  and  moist  adiabatic  lapse  rate. 
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The  upward  long  wave  radiative  flux  (R/t)  is  computed  by  a  formulation  from  Yamada  (27),  i.e., 

Rtf=£(TT£  +  (l-e)RLl  >  (5) 

To  evaluate  equation  5,  one  requires  an  estimate  of  the  effective  ground  or  surface  temperature, 
Tg.  We  achieve  this  through  the  application  of  M-0  similarity  formulations  (10-14)  and  the 
flux-gradient  hypothesis  as  advocated  by  Dyer  (12)  and  Webb  (14). 


The  similarity -based  expressions  for  the  sensible  and  latent  heat  fluxes  are  defined  by  the 
MARIAH  algorithm  (28)  and  are  expressed  as 


H  =  -c„  pu *  0*  =  - 


cp  pk2Vr(Tr-Tg) 

0m  0h  ln  "  d  ln  "  d 


^0 


'■h 


and 


(6) 


L'E 


12  p  kz  Vr  (qr  —qg) 
-L  p  u*  q*  = - — 


z~d  z-d 

0m  0q  ln - ln - 


^0 


'■h 


(7) 


in  which  cp  is  specific  heat;  p  is  air  density;  k  is  von  Karman's  constant;  Vr  is  the  reference  height 
wind  speed;  L'  is  the  latent  heat  of  vaporization;  q,  and  qg  are  the  reference  height  and  effective 
surface  values  for  specific  humidity,  respectively;  u,,  is  the  friction  velocity  (m/s);  61  is  the 
potential  temperature  scaling  constant  (K);  q,,  is  the  specific  humidity  scaling  constant  (kg/kg); 


z  is  height  above  ground  in  meters;  d  is  displacement  height  in  meters  (-70%  of  the  principal 
roughness  element);  zo  is  the  aerodynamic  (surface)  roughness  height  in  meters  (-14%  of  the 
principal  roughness  element);  Zh  is  the  thermal  roughness  height  in  meters  (-13%  of  the 
aerodynamic  roughness);  and  (j)m,  (fa,  and  (J)q  are  the  dimensionless  wind  shear,  dimensionless 
temperature  lapse  rate,  and  dimensionless  humidity  lapse  rate,  respectively.  Here,  the 
dimensionless  term  for  humidity  is  assumed  to  be  equal  to  the  dimensionless  temperature  lapse, 
i.e.,  (j)q  =  (f>h,  even  though  field  observations  have  shown  that  this  is  not  always  correct.  Also,  for 
stable  and  very  stable  conditions  at  night,  Mahrt  (29)  recommends  setting  zo  =  z.i,.  Equations  6 
and  7  can  be  evaluated  in  a  straightforward  manner,  given  expressions  for  the  dimensionless 
shear  and  lapse  rate  terms. 


Following  the  recommendations  of  Hogstrom  (30),  we  use 

0m  =[l-19(z/L)r1/4  and  </>h  =  0.95[l-11.6(z/L)]“1/2  ,  (8) 

for  daytime  unstable  atmospheric  conditions,  and  we  use 

0m  =1  +  5.3 (z/L)  and  (ph  =1  +  8 (z/L)  (9) 
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for  nighttime  stable  atmospheric  conditions',  where  Busch  (11)  defines  the  M-0  scaling  ratio  as 


z_ 

L 


g  Ay* 

ev  u  * 


(10) 


in  which  0V  =  T,(  1  +  0.6 1  qr)  is  the  virtual  potential  temperature  and  #,*  =  0*  +  0.6 1  Trq*  is  the 
virtual  potential  temperature  scaling  constant  ( q *  is  the  specific  humidity  scaling  constant). 
Here,  reference  level  specific  humidity,  qr,  is  determined  as  a  function  of  relative  humidity  (Rh), 
atmospheric  pressure  (Pr),  and  air  temperature  ( Tr)  with  the  formulation  from  Rogers  (31),  i.e., 


3.8  Ru 

Qr  =  e*P 


5.44xl0: 


V 


1 


1 


A 


(11) 


y2.73.15  TrJJ 

In  contrast,  surface  specific  humidity,  qg,  is  determined  as  a  function  of  the  following  parameters 
(32):  reference  level  specific  humidity  (qr),  saturation  specific  humidity  at  the  surface  (qs)  (i.e., 
substitute  into  equation  1 1  Ri,  =  100%  and  T,  =  77),  the  water  capacity  of  the  field  (Of),  and 
finally,  soil  water  content  (wy).  As  a  result,  surface  specific  humidity  can  be  expressed  as 


in  which 


Q g  (l  Qr  "t"  8  Qs  ’ 


(12) 


8  =  0.25 


1-cos 


f  _  AA2 

n  we 


0 
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and  8=  1  for  ws  >  Of  . 

Next,  we  substitute  equations  6  and  7  for  H  and  L  'E  into  equation  1  and  define  the  soil  heat  flux 
(G)  as  G  =  (iRnd  (where  Rnet  =  RlI  -  Rl t  +  Rsi  -  R.sl  and  Rst  =  oRsi).  Then  after  some 
rearranging,  we  solve  equation  1  for  Tg,  which  gives  the  quartic  form  for  the  equilibrium  model 
as 


(1  -/3)(ecr  T*  - (1  - a)Rsl  -  e Ra)+  Tg 


Cp  p  k  Vr 


,  z-d  z-d 
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in  which  a  is  the  surface  short  wave  reflectivity,  commonly  called  “albedo,”  and  /7  i  s  the 
percentage  (e.g.,  10%  daytime  and  5%  nighttime)  of  the  computed  net  radiation  used  to  calculate 
the  soil  (ground  storage)  heat  flux.  Naturally,  at  night,  Rsi  =  0.  The  quartic  equation  is  resolved 
numerically  via  a  Newton-Raphson  iterative  scheme. 


Alternately,  Webb  (14)  suggests  <z>„,  =  $,=  !+  5(z/L )  for  nighttime  stable  conditions. 
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2.2  The  Energy  Budget  Model  Implemented  in  ABF  A 

Wilson  (15)  and  Wilson  and  Szeto  (16)  report  on  the  development  of  the  prototype  model  ABFA 
for  the  analysis  of  atmospheric  refraction  and  turbulence  effects  on  acoustical  detection  and 
tracking  systems.  Among  the  meteorological  input  capabilities  embedded  in  ABFA  is  a  surface 
energy  budget  model  for  generating  low  atmosphere  profiles  of  effective  sound  speed.  The  ABFA 
energy  budget  model  is  based  on  the  empirical  formulations  of  Burridge  and  Gadd  (33),  as 
discussed  in  Stull’s  (21)  textbook  (Chapter  7).  For  the  nighttime  case,  the  net  long  wave  radiation  is 
given  by  the  following  expression: 

/?LT-/?Li  =  (0-08/»pXl-0-0-l<TH -O-Som -0.6<tl)  ,  (15) 

in  which  crH,  (JM,  and  are  the  cloud  cover  fractions  for  high,  middle,  and  low  clouds, 
respectively.  (In  contrast,  cloud  cover  is  accounted  for  in  the  parameterization  for  RlI  shown  in 
equation  4.)  Clearly,  this  is  a  much  simpler  formulation  than  the  solution  of  the  quartic. 

In  ABFA,  the  net  surface  radiation  is  balanced  entirely  by  the  surface  heat  flux,  Q  ,  so  that  the 
M-0  temperature  scaling  constant  is  calculated  as 


6* 


and  the  M-0  scaling  length  is  derived  as 


Q 

pc  pll* 


(16) 


L  =  - 


u 


3 

* 


fc#  JL 

e  pc  p 


(17) 


where  the  friction  velocity  can  be  determined  from  the  reference  height  wind  speed  (Vr)  and 
values  for  roughness  length  (zo)  and  displacement  (d),  i.e., 


u* 


kv, 


fZ~d ' 


(18) 


V  ZO  J 

Thus,  the  ABFA  energy  budget  model  output  is  used  to  estimate  the  similarity-based,  log-linear 
profiles  for  wind  speed  and  temperature  extended  upward  through  the  surface  layer,  as  discussed 
next. 


2.3  The  Similarity-Based  Profiles 

Given  values  for  the  M-0  scaling  constants  from  energy  budget  considerations,  the  formulations 
of  Benoit  (34)  and  Nickerson  and  Smiley  (35)  are  applied  to  determine  the  similarity-based 
profiles  for  wind  velocity,  temperature,  and  specific  humidity  as  follows: 
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where,  for  the  nighttime  case,  the  diabatic  influence  functions  for  momentum  and  heat,  y/m  and 
y/h  respectively,  are  given  by  Webb  (14),  i.e., 


¥„ 


Vh=~  5J  • 


(22) 


Here,  we  assume  that  y/q  =  y/h. 

For  extended  profiles  of  relative  humidity,  we  use  equation  1 1  and  the  following  expression  for 
atmospheric  pressure  given  by  Campbell  and  Norman  (36),  i.e., 
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3.  Model  Results 


3.1  A  Hypothetical  Case 

A  hypothetical  nighttime  case  is  outlined  in  Table  1.  Input  parameter  values  were  chosen  to 
represent  low  wind  speeds  and  cool  temperatures  over  a  relatively  smooth  surface.  For  this 
example,  the  primary  output  parameters  from  the  alternate  energy  budget  model  (described 
previously)  are  given  in  Table  2.  First,  it  is  shown  that  fairly  good  closure  exists  for  this 
calculation,  i.e.,  Rnet  -  G  is  approximately  equal  to  H  +  L  E.  This  result  helps  to  validate  the 

model  for  partitioning  the  net  radiation  into  its  surface  flux  components.  The  small  error  (i.e., 

_2 

A  =  6.83  Wm  “)  may  be  attributed  to  the  approximation  for  soil  heat  flux.  Second,  the  estimate 
for  surface  temperature,  i.e.,  Tg  =  Tr  -  2.25  K,  appears  to  be  quite  reasonable  for  moderately 
stable  conditions  (i.e.,  L  =  35  m).  Third,  the  M-0  scaling  constants  appear  to  be  valid  and  useful, 
so  that  applying  equations  19  through  24  provides  the  nighttime  profiles  for  wind  velocity,  air 
temperature,  and  relative  humidity  for  the  layer  z  <  60  m  (see  Figure  2). 
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Table  1.  Energy  budget  model  parameters 


Parameter 

Symbol 

Value 

Reference  atmospheric  pressure 

Pr 

1000  mbar 

Reference  relative  humidity 

Rh 

40% 

Reference  air  temperature 

Tr 

283.15  K 

Reference  wind  velocity 

Vr 

3.0  ms-1 

Reference  height  above  ground 

Zr 

2.0  m 

Principal  roughness  element  or  canopy 
height 

h 

0.50  m 

Surface  (long  wave)  emissivity 

8 

0.98 

Soil  water  content 

0. 1 3  m3m-3 

Table  2.  Energy  budget  model  results 


Parameter 

Symbol 

Value 

Net  radiative  flux 

R net 

-86.62  Wrn-2 

Soil  (ground  storage)  heat  flux 

G 

-4.33  Wm-2 

Sensible  heat  flux 

H 

-114.50  Wm-2 

Evaporative  (latent)  heat  flux 

L'E 

+  25.38  Wm-2 

Friction  velocity 

u*. 

0.35  ms-1 

Temperature  scaling  constant 

9* 

0.27  °K 

Humidity  scaling  constant 

q* 

-2.35E-05  kgkg-1 

Monin-Obukhov  scaling  length 

L 

34.64  m 

Surface  air  temperature 

Tg 

280.90  °K 

60- 

a) 

60- 

b) 

60- 

c) 

50- 

50- 

50- 

40- 

40- 

40- 

Height  (m) 
8 

_ , _ i _ 

/ 

30- 

30- 

\ 

20- 

/ 

20- 

/ 

20- 

\ 

10- 

J 

10- 

I 

10- 

\ 

0- 

0- 

0- 

2  4  6  10  12  14  30  35  40 

V\lnd  velocity  (ms1)  Air  temperature  (C)  Relative  humidity  (%) 

Figure  2.  Typical  similarity-based  profiles  for  the  nighttime  case  derived  from  the 
energy  budget  model:  a)  wind  velocity  (in  units  ms-1),  b)  air  temperature 
(in  units  C),  and  c)  relative  humidity  (in  percent). 
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Next,  we  apply  equations  1  and  2  and  the  data  shown  before  to  derive  the  extended  profiles  for  the 
effective  speed  of  sound  (see  Figure  3).  In  Figure  3,  the  effective  sound  speed  is  determined  for 
upwind  propagation  (left),  i.e.,  the  bearing  of  the  receiver  from  the  source  and  the  direction  of 
the  wind  from  north  are  the  same,  and  for  downwind  propagation  (right),  i.e.,  the  bearing  angles 
are  180°  opposite.  An  asymmetry  exists  in  the  graph  so  that  the  profile  for  downwind 
propagation  (on  the  right-hand  side)  is  increasing  with  height  above  ground  level,  while  the 
curve  for  upwind  propagation  (on  the  left  side)  is  decreasing  (slightly)  with  height.  The  effect 
that  such  profile  variations  in  sound  speed  will  have  on  short  range  acoustic  attenuation  is 
discussed  next. 


Figure  3.  Low  atmosphere  profiles  of  the  effective  speed  of 
sound  derived  from  the  energy  budget  model  for 
upwind  propagation  (dashed)  and  downwind 
propagation  (solid). 

3.2  Approximation  of  Short  Range  Acoustic  Attenuation 

To  briefly  examine  short  range  acoustic  attenuation  at  night,  we  use  the  low  atmosphere  profiles 
of  wind  speed,  temperature,  and  relative  humidity  (shown  before)  as  input  to  a  flat  earth,  non- 
turbulent  acoustic  propagation  model  called  the  Windows  (version)  Scanning  Fast  Field  Program 
(WSCAFFIP).  WSCAFFIP  is  a  numerical  code  developed  for  assessing  environmental  effects  on 
short  range  acoustic  attenuation  (7,38).  WSCAFFIP  determines  acoustic  attenuation  as  relative 
sound  pressure  loss  with  range  and  azimuth  for  a  given  frequency  and  source-to-receiver  geometry. 
WSCAFFIP  contains  propagation  algorithms  to  represent  the  effects  of  atmospheric  refraction, 
diffraction,  absorption,  and  reflection  (ground  impedance)  on  acoustic  transmission.  Table  3  lists 
the  model  parameters  for  an  initial  approximation  of  short  range  acoustic  attenuation  over  an  open 
grass-covered  ( h  =  0.5  m)  field.  Figures  4  and  5  show  the  WSCAFFIP  results  corresponding  to  the 
modeled  profiles  of  effective  sound  speed  generated  by  the  alternate  (quartic)  model. 
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Table  3.  WSCAFFIP  model  parameters 


Parameter 

Symbol 

Value 

Distance  of  receiver  from  source 

X 

1000  m 

Range  resolution 

Ax 

10m 

Bearing  of  receiver  from  source 

Or 

90° 

Bearing  of  the  wind  from  North 

0W 

90°  (upwind)  and  270°  (downwind) 

Source  height  above  ground 

hs 

2  m 

Receiver  height  above  ground 

h/t 

1  m 

Frequency  of  interest 

f 

100-  1000  Hz 

Ground  porosity  (tall  grass  pasture) 

QP 

0.675 

Flow  resisitivity  (tall  grass  pasture) 

ar 

330  kPa  s  m-2 

Figure  4.  Line  plots  of  acoustic  attenuation  versus  range  for 
upwind  propagation  (dashed)  and  downwind 
propagation  (solid),  (i.e.,  WSCAFFIP  model  results 
for  300  Hz  and  600  Hz,  only). 
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As  stated  earlier,  acoustic  waves  will  tend  to  be  refracted  downward  as  the  effective  sound  speed 
increases  with  height  and  refracted  upward  as  the  effective  sound  speed  decreases  with  height. 
Therefore,  one  might  expect  the  amounts  of  attenuation  (in  part)  to  depend  on  the  strength  of  the 
sound  speed  profile  gradient  between  the  ground  and  the  model  top.  Because  generally,  the 
effective  sound  speed  profile  for  upwind  propagation  has  values  decreasing  with  height,  we  would 
expect  to  see  greater  attenuation  with  range  at  each  frequency  for  upwind  propagation  as  a  result. 
The  line  plots  of  relative  sound  loss  with  range  in  Figure  4  (for  300  Hz  and  600  Hz  only)  and  the 
contour  graphs  of  attenuation  as  a  function  of  range  and  frequency  shown  in  Figure  5  (for  100  to 
1000  Hz  inclusive)  show  this  result  satisfactorily. 

In  contrast,  the  behavior  of  the  downwind  propagation  attenuation  curves  can  be  understood  in 
terms  of  “ducting”  of  acoustic  modes.  For  this  example,  the  sound  speed  profile  forms  an 
acoustic  duct  between  the  surface  and  60  m.  At  lower  frequencies  (100  to  200  Hz),  acoustic 
waves  may  be  too  long  for  any  significant  ducting  to  occur.  At  300  Hz,  there  appear  to  be  one  or 
two  trapped  modes,  while  at  400  Hz  (and  greater)  several  modes  appear  to  be  trapped,  creating 
the  interference  pattern.  In  this  way,  ducting  of  acoustic  energy  across  the  lower  boundary  layer 
at  night  is  greatly  affected  by  local  meteorological  profile  structure. 

3.3  Comparison  to  Observations 

At  the  start  of  this  project,  it  was  thought  feasible  to  present  modeled  profiles  of  effective  sound 
speed  for  the  nighttime  case  in  comparison  to  values  of  effective  speed  of  sound  derived  from 
selected  micro-meteorological  data.  As  an  example,  the  tower  data  used  to  generate  the  three 
time-height  series  shown  in  Figure  6  were  collected  as  part  of  the  1999  Cooperative  Atmosphere- 
Surface  Exchange  Study  (CASES-99)  (39).  CASES-99  was  an  extensive  field  experiment  that 
focused  on  the  nocturnal  boundary  layer  (NBL)  under  mainly  clear  sky,  light  near  surface  wind 
conditions  (this  included  occurrences  of  weak  and  intermittent  turbulence).  The  experimental 
period  was  from  October  1  to  31,  1999.  The  test  area  for  CASES-99  was  near  Leon,  KS 
(50  km)  east  of  Wichita,  KS.  The  test  site  was  chosen  (in  part)  because  of  its  relatively  flat 
terrain  (average  slopes  were  0.5  degree).  A  large  number  of  instruments  were  deployed  in  an 
area  4.8  x  3.2  km.  These  instruments  included  a  heavily  instrumented  55-m  tower,  numerous 
10-m  towers,  multiple  radar,  lidar,  scintillometers,  tethersondes,  rawinsondes,  and  research 
aircraft  (see  http://www.joss.ucar.edu/cases99). 

The  subset  of  CASES-99  data  shown  in  Figure  6  is  from  the  first  intensive  observing  period 
(IOP)  on  5  October  1999  and  includes  the  measurements  of  wind  velocity,  air  temperature, 
relative  humidity  (5  min  averaged)  taken  on  the  55-m  main  tower  approximately  5  km  southeast 
of  Leon,  KS  (latitude:  37  deg  38.88'  N,  longitude:  96  deg  44.14'  W,  elevation  433.7  m  above  sea 
level)  in  an  open  field  of  mostly  dry  grasses,  approximately  0.15  m  in  height  (on  average). 
Surface  roughness  height  for  the  site  was  determined  to  be  about  2  cm  and  thermal  emissivity 
was  estimated  at  0.98  (40). 
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_ Local  time  (hrs) _ | 

Figure  6.  Time-height  series  of  a)  air  temperature  (in  units 
C),  b)  mean  horizontal  wind  speed  (in  units  ms-1), 
and  c)  effective  sound  speed  (in  units  ms-1), 
generated  from  CASES-99  (IOP-1)  field  data. 


However,  we  experienced  some  difficulties  in  the  application  of  traditional  M-0  profile  theory 
for  calculating  nighttime  turbulence  structure  from  temperature,  humidity,  and  wind  speed 
difference  data  (see  equations  6,  7,  and  14).  As  an  example,  Figure  7  shows  a  time  series  of  the 
modeled  surface  temperature  (77,),  the  5-m  measured  temperature  (77),  and  the  5-m  measured 
wind  speed  ( Ur )  for  the  first  nighttime  IOP.  For  this  example,  reference  height  is  zr  =  5  m  and 
the  canopy  height  is  h  =  0.15  m.  Certainly,  the  surface  temperature  estimates  are  unrealistically 
low.  This  “runaway  cooling,”  according  to  Mahrt  (29),  is  a  known  defect  in  computer  models  that 
attempt  to  estimate  meteorological  gradients  near  the  ground  at  night  (during  very  stable  conditions) 
with  M-0  similarity  functions.  To  illustrate  the  problem  further,  Figure  8  shows  a  plot  (at  t  = 

0.125  hr  past  midnight)  of  the  estimated  surface  temperature  (77,),  the  estimated  M-0  scaling  length 
(L),  and  the  estimated  dimensionless  wind  shear  (</),„)  as  a  function  of  the  5-m  reference  height  wind 
speed  (Ur),  as  the  wind  speed  is  varied  from  Ur  =1.0  ms-1  to  Ur  =  10.0  ms1.  In  addition,  the  data  in 
Table  4  and  Figure  9  show  a  more  detail  view  of  the  results  plotted  in  Figure  8.  The  data  clearly 
show  that  the  model  calculation  collapses  at  Ur  =  3.6  -3.8  ms1,  where  Tg  <  5.0  C,  L  <  8.0  m,  and 
(pm  >  1.6.  Then,  for  still  lower  wind  speeds,  L  becomes  very  small — small  enough  to  cause  a 
computational  sensitivity,  as  evidenced  in  the  plot  for  (pm.  Given  that  both  (pm  and  (ph  are  needed  to 
define  two  principal  components  of  the  quartic,  it  can  be  reasoned  that  when  L  is  small  and  the  (p' s 
are  large,  the  model  defaults  to  a  kind  of  radiative  equilibrium,  where  little  heat  energy  (sensible, 
latent,  or  released  from  the  soil)  is  available  to  offset  the  radiative  term.  When  this  occurs,  i.e., 
when  turbulence  collapses,  the  model  results  are  invalid. 

As  a  final  point,  we  consider  the  effect  on  the  model  results  brought  about  by  varied  surface 
roughness  (see  Figure  10).  Recall  that  for  this  example,  the  canopy  height  is  h=0.15  m  (short 
grass)  and  the  roughness  length  is  zo  =  0.02  m.  In  Figure  10,  when  Ur  =  3.6  ms1  and  the  canopy 
height  <  0.25  m,  the  model  again  breaks  down.  Thus,  for  whatever  the  order  of  wind  speed  and 
roughness  length,  there  appears  be  a  derivable  limit  beyond  which  the  calculation  should  be 
discounted.  The  arrow  on  Figure  9  suggests  such  a  limit. 

It  is  interesting  to  note  that  others  had  previously  reported  difficulties  in  the  use  of  similarity- 
based  models  for  the  nighttime  case  (in  the  calculation  of  optical  turbulence  intensity  [  \ ) 

(41-43).  As  an  example,  Figure  11  shows  results  from  the  CN2  over-land  optical  turbulence 
model  (42,43),  derived  from  two  levels  (above  ground)  of  conventional  micro-meteorological 
information  as  input  (i.e.,  wind  speed,  temperature,  and  humidity  retrieved  on  farmland,  over  flat 
barren  soil,  under  mostly  clear  skies).  Although  the  CN2  model  estimates  are  in  fairly  good 
agreement  with  the  observations  during  daytime  hours,  there  are  several  instances  at  night  when 
the  model  results  show  either  a  consistent  bias  (about  a  factor  of  2  and  mainly  underestimated)  in 
comparison  to  the  measured  turbulence  data  or  the  model  results  fluctuate  toward  extreme 
(minimum)  values.  Here,  the  largest  errors  occur  when  there  is  a  computational  sensitivity,  as 
described  before.  Average  surface  layer  wind  speeds  for  the  nighttime  periods  shown  in 
Figure  11  were  approximately  3  to  4  ms1. 
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Figure  7.  Time  series  of  a)  modeled  surface  temperature  ( Tg )  and 
the  5-m  temperature  (Tr)  and  b)  the  5-m  measured 
wind  speed  (U,)  for  the  first  nighttime  IOP.  (The  dotted 
line  in  Figure  7b  helps  to  show  the  large  amount  of  data 
in  this  example  where  U  <  2.0  ms-1.) 
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Figure  8.  A  plot  at  t  =  0.125  hr  of  a)  the  estimated  surface  temperature 
( Tg ),  b)  the  estimated  M-O  scaling  length  (L),  and  c)  the 
estimated  dimensionless  wind  shear  (</>„,)  as  a  function  of  the 
5-m  reference  height  wind  speed  ( Ur ). 


Table  4.  Breakdown  of  M-O  similarity  models  for  the  nighttime  case 


U r  (ms  ') 

Tg  <c) 

L  (m) 

3.2 

-8.698 

0.005 

837.346 

3.3 

-8.698 

0.008 

547.225 

3.4 

-8.697 

0.016 

275.226 

3.5 

-8.682 

0.066 

71.424 

3.6 

-8.447 

0.275 

18.173 

3.7 

2.579 

4.914 

1.964 

3.8 

5.140 

8.477 

1.559 

3.9 

6.062 

10.945 

1.433 

4.0 

6.647 

13.245 

1.358 

4.1 

7.072 

15.504 

1.306 

4.2 

7.403 

17.780 

1.267 

4.3 

7.671 

20.096 

1.236 

4.4 

7.896 

22.471 

1.211 

4.5 

8.088 

24.908 

1.190 

4.6 

8.255 

27.426 

1.173 

4.7 

8.401 

30.016 

1.158 

4.8 

8.532 

32.693 

1.145 
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Figure  1 1 .  Time  series  of  CN2  model  output  (dashed  circle)  compared  to 

scintillometer  data  (solid  square),  determined  from  input  data  at  1 
and  4  m  above  ground  collected  over  a  450-m  path  on  a)  8-9  July 
1992  and  b)  13-14  July  1992  (42,43). 
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4.  Summary  and  Conclusions 


The  Army  is  looking  to  implement  acoustic  propagation  models  that  will  incorporate  the  impacts 
of  the  atmosphere,  turbulence,  and  terrain  on  acoustic  signatures.  A  key  element  in  estimating 
outdoor  acoustic  attenuation  with  range  (and  frequency)  is  modeling  the  variation  of  the  effective 
speed  of  sound  through  the  lower  atmosphere.  M-0  similarity  models  have  been  very  well 
suited  for  this  application,  especially  for  daytime  convective  atmospheres,  e.g.,  the  energy 
budget  and  meteorological  models  in  ABFA.  However,  for  the  nighttime  case,  we  were 
interested  in  testing  the  effectiveness  of  an  earlier  quartic  solution  for  the  surface  temperature, 
i.e.,  an  alternate  energy  budget  model,  which  explicitly  computes  the  net  surface  radiation  and 
the  fluxes  of  heat  and  moisture  close  to  the  ground. 

For  a  hypothetical  nighttime  case  representing  moderately  stable  atmospheric  conditions,  the 
alternate  (quartic)  model  performed  quite  well.  Model  estimates  for  effective  surface 
temperature,  net  radiation,  sensible  heat  flux,  soil  heat  flux,  and  latent  heat  flux  were  comparable 
to  earlier  published  results  (17,18).  Also,  for  this  case,  energy  budget  closure  was  met 
satisfactorily.  Thus,  the  primary  energy  budget  model  results  could  be  used  to  predict  profiles  for 
wind  speed,  air  temperature,  and  relative  humidity,  from  which,  effective  sound  speed  (and 
ultimately  acoustic  attenuation)  information  could  be  generated. 

However,  we  found  that  reasonable  values  for  effective  surface  temperature  were  not  obtainable 
from  the  alternate  model  to  simulate  the  first  nighttime  IOP  of  the  CASES-99  study. 
Computational  difficulties  were  experienced  because  of  low  roughness  height  and  low  wind 
speed.  We  found  that  as  the  M-0  scaling  length  becomes  small,  the  energy  budget  model  breaks 
down  and  the  computation  is  not  viable.  Thus,  until  better  solutions  can  be  realized  through 
future  modeling  or  experimental  works,  we  recommend  that  the  alternate  model  be  applied  only 
when  input  canopy  height  and  reference  level  wind  speeds  can  affect  values  for  the  M-0  scaling 
length  larger  than  about  13  m. 
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